Evaluation of SMAP Freeze/Thaw Retrieval Accuracy at Core Validation Sites in the Contiguous United States

Seasonal freeze-thaw (FT) impacts much of the northern hemisphere and is an important control on its water, energy, and carbon cycle. Although FT in natural environments extends south of 45°N, FT studies using the L-band have so far been restricted to boreal or greater latitudes. This study addresses this gap by applying a seasonal threshold algorithm to Soil Moisture Active Passive (SMAP) data (L3_SM_P) to obtain a FT product south of 45°N (‘SMAP FT’), which is then evaluated at SMAP core validation sites (CVS) located in the contiguous United States (CONUS). SMAP landscape FT retrievals are usually in good agreement with 0–5 cm soil temperature at SMAP grids containing CVS stations (\u3e70%). The accuracy could be further improved by taking into account specific overpass time (PM), the grid-specific seasonal scaling factor, the data aggregation method, and the sampling error. Annual SMAP FT extent maps compared to modeled soil temperatures derived from the Goddard Earth Observing System Model Version 5 (GEOS-5) show that seasonal FT in CONUS extends to latitudes of about 35–40°N, and that FT varies substantially in space and by year. In general, spatial and temporal trends between SMAP and modeled FT were similar

Application of observation operators for field scale soil moisture averages and variances in agricultural landscapes

Scale difference between in situ and remotely sensed soil moisture observations and model grid size has been an issue for validation of remote sensing data, soil moisture data assimilation and calibration of hydrologic models. This study aims to link two different scales of soil moisture estimates by upscaling single point measurements to field averages for representing field-scale agricultural areas (∼2 ha) located within the Upper Cedar Creek Watershed in northeastern Indiana. Several statistical methods, mainly focusing on cumulative distribution function (CDF) matching, are tested to upscale point measurements to spatially representative soil moisture. These transforming equations are termed observation operators. The CDF matching is found to be the most successful upscaling method followed by the mean relative difference method using temporally stable point measurements. This study also tests the temporal and spatial (horizontal and vertical) transferability of the different observation operators. Results indicate that the two observation operators from the CDF matching approach and the mean relative difference method using a temporally stable location are transferable in space, but not in time. Rainfall characteristic is most likely the dominant factor affecting the success of observation operator transferability. In addition, the CDF matching approach is found to be an effective method to deduce spatial variability (standard deviation) of soil moisture from single point measurements

Combined Radar-Radiometer Surface Soil Moisture and Roughness Estimation

A robust physics-based combined radar-radiometer, or Active-Passive, surface soil moisture and roughness estimation methodology is presented. Soil moisture and roughness retrieval is performed via optimization, i.e., minimization, of a joint objective function which constrains similar resolution radar and radiometer observations simultaneously. A data-driven and noise-dependent regularization term has also been developed to automatically regularize and balance corresponding radar and radiometer contributions to achieve optimal soil moisture retrievals. It is shown that in order to compensate for measurement and observation noise, as well as forward model inaccuracies, in combined radar-radiometer estimation surface roughness can be considered a free parameter. Extensive Monte-Carlo numerical simulations and assessment using field data have been performed to both evaluate the algorithms performance and to demonstrate soil moisture estimation. Unbiased root mean squared errors (RMSE) range from 0.18 to 0.03 cm3cm3 for two different land cover types of corn and soybean. In summary, in the context of soil moisture retrieval, the importance of consistent forward emission and scattering development is discussed and presented

Characterization of Forest Opacity Using Multi-Angular Emission and Backscatter Data

This paper discusses the results from a series of field experiments using ground-based L-band microwave active/passive sensors. Three independent approaches are employed to the microwave data to determine vegetation opacity of coniferous trees. First, a zero-order radiative transfer model is fitted to multi-angular microwave emissivity data in a least-square sense to provide "effective" vegetation optical depth. Second, a ratio between radar backscatter measurements with the corner reflector under trees and in an open area is calculated to obtain "measured" tree propagation characteristics. Finally, the "theoretical" propagation constant is determined by forward scattering theorem using detailed measurements of size/angle distributions and dielectric constants of the tree constituents (trunk, branches, and needles). The results indicate that "effective" values underestimate attenuation values compared to both "theoretical" and "measured" values

Physics-Based Retrieval of Surface Roughness Parameters for Bare Soils from Combined Active-Passive Microwave Signatures

In the past the effect of soil roughness was often considered secondary within the determination of soil moisture from remote sensing data. Several studies showed that accurate determination of soil roughness leads to an improved estimation of soil moisture. Two default parameters to describe the surface roughness are the standard deviation of the surface height variation and the surface correlation length with its corresponding autocorrelation function. Both parameters (,) affect the emissivity measured by radiometers as well as the backscattering observed by radars. In this study, we develop a physics-based approach to retrieve and by combining both microwave signals based on active-passive microwave covariation. To test the approach, containing a forward model and a retrieval algorithm, we used active/passive microwave data measured with the ComRAD truck-based SMAP simulator at L-band. Results and validations with corresponding field measurements on ground show that and can be estimated simultaneously when using this approach. The physics-based retrieval algorithm works robustly for two investigated test fields having an RMS-Error of 0.68 cm and 0.69 cm between the microwave-based and field-measured -values, and of 3.13 cm and 3.04 cm for -values. The first validation of the results reveals that the influence of the autocorrelation function, needed within the retrieval, is distinct

Impact of Conifer Forest Litter on Microwave Emission at L-Band

This study reports on the utilization of microwave modeling, together with ground truth, and L-band (1.4-GHz) brightness temperatures to investigate the passive microwave characteristics of a conifer forest floor. The microwave data were acquired over a natural Virginia Pine forest in Maryland by a ground-based microwave active/passive instrument system in 2008/2009. Ground measurements of the tree biophysical parameters and forest floor characteristics were obtained during the field campaign. The test site consisted of medium-sized evergreen conifers with an average height of 12 m and average diameters at breast height of 12.6 cm. The site is a typical pine forest site in that there is a surface layer of loose debris/needles and an organic transition layer above the mineral soil. In an effort to characterize and model the impact of the surface litter layer, an experiment was conducted on a day with wet soil conditions, which involved removal of the surface litter layer from one half of the test site while keeping the other half undisturbed. The observations showed detectable decrease in emissivity for both polarizations after the surface litter layer was removed. A first-order radiative transfer model of the forest stands including the multilayer nature of the forest floor in conjunction with the ground truth data are used to compute forest emission. The model calculations reproduced the major features of the experimental data over the entire duration, which included the effects of surface litter and ground moisture content on overall emission. Both theory and experimental results confirm that the litter layer increases the observed canopy brightness temperature and obscure the soil emission

Effective Tree Scattering at L-Band

For routine microwave Soil Moisture (SM) retrieval through vegetation, the tau-omega [1] model [zero-order Radiative Transfer (RT) solution] is attractive due to its simplicity and eases of inversion and implementation. It is the model used in baseline retrieval algorithms for several planned microwave space missions, such as ESA's Soil Moisture Ocean Salinity (SMOS) mission (launched November 2009) and NASA's Soil Moisture Active Passive (SMAP) mission (to be launched 2014/2015) [2 and 3]. These approaches are adapted for vegetated landscapes with effective vegetation parameters tau and omega by fitting experimental data or simulation outputs of a multiple scattering model [4-7]. The model has been validated over grasslands, agricultural crops, and generally light to moderate vegetation. As the density of vegetation increases, sensitivity to the underlying SM begins to degrade significantly and errors in the retrieved SM increase accordingly. The zero-order model also loses its validity when dense vegetation (i.e. forest, mature corn, etc.) includes scatterers, such as branches and trunks (or stalks in the case of corn), which are large with respect to the wavelength. The tau-omega model (when applied over moderately to densely vegetated landscapes) will need modification (in terms of form or effective parameterization) to enable accurate characterization of vegetation parameters with respect to specific tree types, anisotropic canopy structure, presence of leaves and/or understory. More scattering terms (at least up to first-order at L-band) should be included in the RT solutions for forest canopies [8]. Although not really suitable to forests, a zero-order tau-omega model might be applied to such vegetation canopies with large scatterers, but that equivalent or effective parameters would have to be used [4]. This requires that the effective values (vegetation opacity and single scattering albedo) need to be evaluated (compared) with theoretical definitions of these parameters. In a recent study [9], effective vegetation opacity of coniferous trees was compared with two independent estimates of the same parameter. First, a zero-order RT model was fitted to multiangular microwave emissivity data in a least-square sense to provide effective vegetation optical depth as done in spaceborne retrieval algorithms. Second, a ratio between radar backscatter measurements with a corner reflector under trees and in an open area was calculated to obtain measured tree propagation characteristics. Finally, the theoretical propagation constant was determined by forward scattering theorem using detailed measurements of size/angle distributions and dielectric constants of the tree constituents (trunk, branches, and needles). Results indicated that the effective attenuation values are smaller than but of similar magnitude to both the theoretical and measured values. This study will complement the previous work [9] and will focus on characterization of effective scattering albedo by assuming that effective vegetation opacity is same as theoretical opacity. The resultant effective albedo will not be the albedo of single forest canopy element anymore, but it becomes a global parameter, which depends on all the processes taking place within the canopy including multiple scattering as described

High-Resolution Soil-Moisture Maps Over Landslide Regions in Northern California Grassland Derived From SAR Backscattering Coefficients

Slow-moving landslides are destabilized by accumulated precipitation and consequent soil moisture. Yet, the continuous high-resolution soil-moisture measurements needed to aid the understanding of landslide processes are generally absent in steep terrain. Here, we produce soil-moisture time-series maps for a seasonally active grassland landslide in the northern California coast ranges, USA, using backscattering coefficients from NASA's uninhabited aerial vehicle synthetic aperture radar at 6-m resolution. A physically based radar scattering model is used to retrieve the near-surface (5-cm depth) soil moisture for the landslide. Both forward modeling (backscattering estimation) and the retrieval (soil-moisture validation) show good agreement. The root-mean-square errors (RMSE) for vertical transmit vertical receive (VV) and horizontal transmit horizontal receive (HH) polarizations in forward model comparison are 1.93 dB and 1.88 dB, respectively. The soil-moisture retrieval shows unbiased RMSE of 0.054 m³/m³. Our successful retrieval benefits from the surface and double-bounce scattering, which is common in grasslands. The retrieved maps show saturated wetness conditions within the active landslide boundaries. We also performed sensitivity tests for incidence angle and found that the retrieval is weakly dependent on the angle, especially while using copolarized HH and VV together. Using the two copolarized inputs, the retrieval is also not sensitive to the change of orientation angles of grass cylinders. The physical model inversion presented here can be generally applied for soil-moisture retrieval in areas with the same vegetation cover types in California

Comparison between Dense L-Band and C-Band Synthetic Aperture Radar (SAR) Time Series for Crop Area Mapping over a NISAR Calibration-Validation Site

Crop area mapping is important for tracking agricultural production and supporting food security. Spaceborne approaches using synthetic aperture radar (SAR) now allow for mapping crop area at moderate spatial and temporal resolutions. Multi-frequency SAR data is highly useful for crop monitoring because backscatter response from vegetation canopies is wavelength dependent. This study evaluates the utility of C-band Sentinel-1B (Sentinel-1) and L-band ALOS-2 (PALSAR) data, collected during the 2019 growing season, for generating accurate active crop extent (crop vs. non-crop) classifications over an agricultural region in western Canada. Evaluations were performed against the Agriculture and Agri-Food Canada satellite-based Annual Cropland Inventory (ACI), an open data product that maps land cover across the extent of Canada’s agricultural land. Classifications were performed using the temporal coefficient of variation (CV) approach, where an optimal crop/non-crop delineating CV threshold (CVthr) is selected according to Youden’s J-statistic. Results show that crop area mapping agreed better with the ACI when using Sentinel-1 data (83.5%) compared to PALSAR (73.2%). Analysis of performance by crop reveals that PALSAR’s poorer performance can be attributed to soybean, urban, grassland, and pasture ACI classes. This study also compared CV values to in situ wet biomass data for canola and soybeans, showing that crops with lower biomass (soybean) had correspondingly lower CV value

Partitioning Evapotranspiration in Semiarid Grassland and Shrubland Ecosystems Using Diurnal Surface Temperature Variation

The encroachment of woody plants in grasslands across the Western U.S. will affect soil water availability by altering the contributions of evaporation (E) and transpiration (T) to total evapotranspiration (ET). To study this phenomenon, a network of flux stations is in place to measure ET in grass- and shrub-dominated ecosystems throughout the Western U.S. A method is described and tested here to partition the daily measurements of ET into E and T based on diurnal surface temperature variations of the soil and standard energy balance theory. The difference between the mid-afternoon and pre-dawn soil surface temperature, termed Apparent Thermal Inertia (I(sub A)), was used to identify days when E was negligible, and thus, ET=T. For other days, a three-step procedure based on energy balance equations was used to estimate Qe contributions of daily E and T to total daily ET. The method was tested at Walnut Gulch Experimental Watershed in southeast Arizona based on Bowen ratio estimates of ET and continuous measurements of surface temperature with an infrared thermometer (IRT) from 2004- 2005, and a second dataset of Bowen ratio, IRT and stem-flow gage measurements in 2003. Results showed that reasonable estimates of daily T were obtained for a multi-year period with ease of operation and minimal cost. With known season-long daily T, E and ET, it is possible to determine the soil water availability associated with grass- and shrub-dominated sites and better understand the hydrologic impact of regional woody plant encroachment